Adiabatic Charge Pumping in Almost Open Dots 
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We consider adiabatic charge transport through an almost open quantum dot. We show that 
the charge transmitted in one cycle is quantized in the limit of vanishing temperature and one- 
electron mean level spacing in the dot. The explicit analytic expression for the pumped charge 
at finite temperature is obtained for spinless electrons. The pumped charge is produced by both 
non-dissipative and dissipative currents. The latter are responsible for the corrections to charge 
quantization which are expressed through the conductance of the system. 

PACS numbers; 73.23.Hk, 73.40. Ei, 72.10.Bg 



Adiabatic charge pumping occurs in a system sub- 
jected to a slow periodic perturbation. Upon the comple- 
tion of the cycle, the Hamiltonian of the system returns 
to its initial form, however, a finite charge may be trans- 
mitted through some cross-section of the system. The 
natural question is what is the value of this charge trans- 
mitted through the system during one cycle, Q. This 
question does not have a universal answer. Thouless |] 
showed that for certain one dimensional systems with a 
gap in the excitation spectrum in the thermodynamic 
limit the charge Q is quantized. Such quantized charge 
pumping could be of practical importance as a standard 
of electric current . The accuracy of charge quantiza- 
tion depends on the degree of adiabaticity of the process. 

The practical attempts at creating a quantized electron 
pump use a different approach based on the phenomenon 
of Coulomb blockade In this kind of devices, one 

uses several single electron transistors (SET) connected 
in series to increase the accuracy of charge quantization. 
At least two SET's are necessary to obtain a non-zero 
charge transfer during one cycle. 

Recently, another family of the Coulomb blockade de- 
vices was demonstrated Q - semiconductor based quan- 
tum dots. The advantage of these devices is the pos- 
sibility of changing not only the gate voltage and thus 
the average electron number in the dot but also the con- 
ductance of the quantum point contacts (QPC's) sep- 
arating the dot from the leads. By doing so, one can 
traverse from almost classical Coulomb blockade to the 
completely open dot where the effects of the charge quan- 
tization are diminished. 

Having those semiconductor structures in mind, we 
theoretically study in this Letter the following setup for 
the adiabatic quantum pump. The device is depicted in 
Fig. 1 and consists of a quantum dot connected to two 
reservoirs labeled by a = ±1 by one channel QPC's char- 
acterized by reflection amplitudes ra{t) and capacitively 
coupled to the metallic gate G. The gate potential is char- 
acterized by the number of electrons N which minimizes 
the electrostatic energy of the dot. Experiments with a 
similar setup were reported in Ref. |^. 




FIG. 1. Schematic drawing of a quantum dot electrostat- 
ically defined on a surface of a two dimensional electron gas. 
The dot is connected to two leads by single channel QPC's 
labeled by a. The voltages on the gates "G" and "±1" de- 
termine respectively the average electron number in the dot, 
N{t), and the refiection amplitudes, r±i(t), in the QPC's. 

We show below that the charge adiabatically transfered 
in one cycle is quantized in the limit of zero temperature 
and zero mean level spacing A. For spinless electrons it 
is given by 

^ = z(t) = ^r.(Oe-^'W. (1) 
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We assume that during the cycle the system does not go 
through the degeneracy point z = 0. The degree of adia- 
baticity of the process depends on the proximity to this 
degeneracy point. 

To illustrate this quantization qualitatively, let us first 
consider the following trivial limit of the pumping cycle: 

i) the right contact (a — -Kl) is completely pinched off, 

ii) the gate voltage is changed from to A^'o, iii) the left 
contact (a — —1) is adiabatically closed. In this pro- 
cess the average charge on the dot changes from A'q to 
the nearest integer rig (If A'q is not a degeneracy point, 
which is assumed, rip is unique), iv) the right contact is 
opened, and the gate voltage is changed from A'o back to 
0, v) the left contact is opened. As a result of this cycle 
the total charge transfered from left to right is an integer 
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no- The non-trivial statement is that even for contacts 
which stay weakly reflecting at all times during the cycle 
the total transfered charge is still quantized and is given 
by Eq. (P. 

We assume that the electrons in the leads are non- 
interacting, and the interaction between the electrons on 
the dot is described by the standard model of Coulomb 
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blockade for strong tunneHng see Sec. 4.3 of Ref. ||] 
for the discussion of the applicability of this model. If the 
dot is connected to the leads of by single channel QPC's 
the system can be treated within the one-dimensional 
effective Hamiltonian |^,^ (we put h = ks = every- 
where) 

Mx)^RAax) ■■ +N{t)\ . (2) 



The first two terms in Eq. (||) describe the non- 
interacting electrons in the left {a = —1) and the right 
(a = +1) QPC respectively: the first term is the lin- 
earized version of the kinetic energy [vp is the Fermi 
velocity) and the second term corresponds to backscat- 
tering in the QPC's. We restrict ourselves to the case 
of small refiection amplitudes, ra{t) ^ 1. Finally, the 
last term describes the effect of the Coulomb blockade in 
the dot, and Ec is its charging energy. For simplicity, 
we explicitly consider the case of spinless electrons, the 
results for the electrons with spins will be presented at 
the end of the paper. The only difference of our Hamilto- 
nian from those considered before resides in its time 
dependence. Since upon the completion of the cycle the 
total charge of the dot returns to its original value the 
integrated current can be calculated through any cross- 
section of the system. For convenience we take half the 
sum of the currents fiowing through the left and right 
point contacts: 



The bosonic operators 4>±^{x) satisfy the following com- 
mutation relations 



cf>Ux),4>Uy) = $l{x),^l{y) 
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The last of Eqs. (g) ensures the correct anticommutation 
relation between left and right moving fermions, how- 
ever, this subtlety will not be important for the problem 
in hand. 

The odd modes 0'l are decoupled from the rest of the 
Hamiltonian and do not contribute to the current (|^), 
hence we can omit them. The relevant part of the Hamil- 
tonian (0) acquires the following form 
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where e is the electron charge. 

Similarly to Ref. the Hamiltonian (|2|) can be 
bosonized according the rules 



where the complex function of time z(t) is given by 
Eq. m). The bosonized current operator (H) becomes 
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where fja are Majorana fermions {fja = ry^, {flajfja'} — 
25a,Q')j whereas the scale A characterizes the large mo- 
mentum cut-off and is of the order of the Fermi wave- 
length. 

Instead of the left and the right modes in Eq. (^) it is 
convenient to introduce the even and odd modes ( 



As one can see from Eq. (|^), the mode 4>'^{0) is pinned 
by the charging energy Ec to the value ttN. Since Eq is 
the large scale in the problem, for the description of the 
low energy dynamics of the system we can integrate (p'^ 
out 1^] and obtain the Hamiltonian of the form 
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H'{t) = 



An 



-oo 



dx 



27r3A 



z{t) exp(i0^(O)) + h.c. 
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where ln7 = C « 0.5772 is the Euler constant. 

The Hamiltonian (|^) reduces to a non-interacting form 
through the introduction of the fermion fields M 
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: exp 
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where C is a Majorana fermion. 

Tire re-fermionized Hamiltonian (^) can be conve- 
niently written in a matrix form 



dxf\x)Hof{x) 
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where k = y'jvpE^jTr^. The re-fermionized current op- 
erator (H) acquires the form 



ip = evp : ^'t(0)*(0) 
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Next, we define the matrix Green function 
G^(t,t';x,x') which includes both normal and anoma- 
lous components 



G{t,t';x,x') = -i(Ttf{x,t)(E)V{x',t')^ 



(13) 



where (. . .) denotes averaging over the quantum state of 
the system. The retarded and advanced Green functions 
G^/"^(i, t'] y, y') are defined in a similar manner. 

All of the observables can be expressed through the 
Green function at almost coinciding times, Q{t]x^x') 



git;x,x') = G{t,t + 0;x,x'), 



(14) 



where G{t,t';x,x') is defined in Eq. (|T^. For example, 
the instantaneous current through the system is given by 



I{t) = -ievpgii{t;x = 0,x' = 0). 



(15) 



In the leading adiabatic approximation the Green func- 
tions of the system coincide with those in equilibrium at 
the instantaneous value of z, and the current (^5|) van- 
ishes. Therefore, to find the current flowing through the 
system in response to an adiabatic change of the Hamilto- 
nian we have to find the first non-adiabatic correction to 
the Green function (H). It obeys the evolution equation 



dg{t) 
dt 
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where i?o(i) is given by Eq. (11c). From this equation 
it follows that the first non-adiabatic correction to the 
Green function ( p^ ) is given by 



sg{t) 



duj 
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dt 
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where ^o[-2(i)] is the equilibrium Green function (|TJ) 
for the instantaneous value of the parameter z{t), and 
Gq^^ [z{t), oj] are the frequency representations for the re- 
tarded and advanced Green functions at fixed z{t). They 
can be found from the following equations 

(co±diag [1, 1, 1/2] - iJo(x)) G^'^{z{t), u; x, x') = 

dieig[S{x-x'),6{x-x'),l], (18a) 
Goizit)]^ f ^npiu;)[G^iz,cu)~G^{z,Lu)), (18b) 



where np{uj) = [1 -I- exp(w/r)] is the Fermi distribu- 
tion function, and lu± = w ± iO. 

We solve Eq. ( |18a|) fo r the Green functions, substitute 
the result in Eqs. (18t) and (p^), and thus find the cur- 
rent (|l^) . Integrating the result over the cycle period we 
find the charge transmitted during the cycle. It can be 
represented through the the dimensionless conductance 
of the system g (in units of /2nT%) as 
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where the dimensionless conductance is given by 
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with C (x, y) being the Riemann zeta-function. At low 
temperatures, T ^ \z\'^Ec, conductance vanishes as 
g (X T'^ / {\z\'^EcY , which means that the transmitted 
charge tends to its quantized value. At high tempera- 
tures T ^ \z\'^Ec the conductance approaches the classi- 
cal value g — 1/2, and the pumped charge ( p^ vanishes. 
In general, for g 0, the pumped charge depends on the 
shape of the contour z{t) and is not a topological number. 

We stress that Eq. (jl^) for the transmitted charge is 
given by the sum of two terms: i) the first one arises from 
non-dissipative currents (This contribution is quantized 
and represents a topological invariant of the cycle.), and 
ii) the second one, containing the conductance, is due to 
dissipative currents generated by the cycling of the dot. 

This fact is not accidental and becomes more trans- 
parent from the following consideration which applies to 
a more general class of systems. The time dependent 
Hamiltonian can be written as H{t) = {t)Ho{t)A{t), 
where Hq is diagonal. In the adiabatic limit the transmit- 
ted charge can be most conveniently evaluated by going 
to the adiabatically rotating basis |V'i(i)) = ^(OlV'j(i)) 
(the "rotating axis representation" ||l^), and calculat- 
ing the current in response to the arising perturbation 
idtA{t)A^ {t) using the Kubo formula 
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dt{[dtAit)A\t)JiO)]). 



(20) 



Now let us apply Eq. ( po| ) to the case of an open dot 
which is connected to the leads by two groups of chan- 
nels denoted by index — ±1, and which is described 
by the Hamiltonian similar to Eq. (j2|). For this purpose 
we introduce the partial particle number operator in each 
channel 



dxipl{aix)'il;i{aix) 



(21) 



We notice, that even though the average particle num- 
ber {fii) H in each channel is infinite, its change during 
the pumping cycle is a well defined quantity and is de- 
termined by the the gate voltage, N{t) and the reflection 
coefficients Vj. The calculation is facilitated by the ex- 
plicit form of the unitary operator 



i(t) =[|cxp 



nTaj{n.j)H 

Vp 



+ 00 



dx Ii{x) 



(22) 



where Ij{x) is the partial particle current operator in 
channel j at point x, 

= ^ [v-j^,^, - V'j] • (23) 

Recalling that the time evolution of operators at a; 7^ 
corresponds to free propagation with velocity vp we read- 
ily express the pumping current through the dimension- 
less partial conductances gtj between the channels 
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H 



dt 



{Sij gij), 
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dit [/,(t,0); 7,(0,0)] . (24) 



For the model of spinless electrons, Eq. (|24| ) repro- 
duces Eq. ([l9|), since is the two terminal con- 
ductance of the system, gn — gi-i = g because of 
charge conservation, and {fii -\- n_i)jj — —N{t), be- 
cause of the large charging energy. The explicit relation 
{fii — h-i)H = — ^lmlnz(i) with z given by Eq. (|l|) re- 
quires the model assumption 

Equation (^) holds even in the case when the elas- 
tic cotunneling ||Tl[| can not be neglected. As shown in 
Ref. joj the elastic returns in the strong tunneling case 
can be described by an iV-independent action. Therefore, 
the explicit form of the unitary operator (^) and Eq. ( |2^ ) 
are still valid. This observation gives us an estimate of 
A/{\z\'^Ec) for the sample-specific correction to 

the quantized value of the charge due to the finite level 
spacing in the dot, A, since the conductance in the val- 
ley of the CB remains finite even in the zero-temperature 
limit g ~ A/(|z|2i?c7) |]. 



Let us now apply Eq. ( p^ ) to strong inelastic cotunnel- 
ing of electrons with spin. Because of the spin symmetry 
and charge conservation we have X)s5q q' ~ 3/2; where 
s is the spin index and g is the two terminal conductance, 
and we readily obtain 



d{ni] 



H 



dt 



[l-5(t)]. 



(25) 



At high temperature t/ — > 1, and pumping is suppressed, 
whereas at low temperatures the conductance vanishes 
as [||. The explicit relation between {n^)H and 
the parameters of the Hamiltonian is a difficult prob- 
lem which was not solved. However, for the anisotropic 
case, where the reflection coefficient in one contact is 
much stronger than in the other, one can show that 
Ea,.arf("a)ffM « (d7V(i)M)sgn(^^a|r„|). More- 
over, if in the transition region, ri ^ r„i, we do not 
go through a degeneracy point costtTV = the pumped 
charge remains quantized at zero temperature, and the 
problem is topologically equivalent to the spinless case. 
All estimates of the finite level spacing effect remain the 
same as for the spinless case. 

To summarize, we considered adiabatic charge trans- 
port through a quantum dot. We have shown, that even 
in the case of small backscattering in the channel, when 
quantum fluctuations of the dot charge are large, the 
transmitted charge is still quantized. We have calcu- 
lated corrections to the quantized value of the charge 
due to finite temperature and the level spacing in the 
dot. These corrections are expressed through the con- 
ductance g and originate from the dissipative currents 
generated by pumping. 
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